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ABSTRACT 


Thermophoresis, the drift of sub-micr on-sized parti- 
cles in tne direction of decreasing temperature, is of interesl 
in many flows with engineering applications It occurs in situa- 
tions involving high temperatures and temperature gradients 
Fe^v researchers have considered the effects of these parameterj 
on the tnermophysical properties of the basic flow This thesis 
incorporates tnese effects in studying the two-dimensiona 
subsonic flow over a cola circular cylinder The approach is 
numerical in nature . finite differenced equations are usee 
to stucy the effects of ^sa^’lous parameters on the flow anc 
on the thermophoresis It is found that the predicted velocity 
temperature, ana concentration profiles are significantly chan- 
ged by the incorporation of variable thermophysical properties 
in the analysis 
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CHAPTER 1 


INTRODUCTION 


Thermophoresxs is the phenomenon by which submicron - sized 
aerosol particles suspended in non-isothermal gas streams acq- 
uire a velocity in the direction ooposite to the temperature 
gradient This phenomenon causes the particles to migrate from 
hotter to colder regions of the gas 

The explanation of this phenomenon is based upon kinetic 
theory arguments, i e , tie particle receives a greater number 
of molecular impacts on the hotter siae than the colder one, 
thus leading to net rate of change of momentum in the direction 
opposite to the temperature gradient The resulting radiometric 
force experienced by the particles is called the therraophoretic 
force. This force, like other interactions between particles 
ana gas, depends essentially upon (1) the Knudsen number, 

(= the ratio of the mean free path, A of the gas molecules to 
the particle radius, a and (2) on the ratio of the thermal 
conductivities of the particles and of surrounding gas Thermo- 
phoresis IS of considerable importance for particles having 
radii as large as 10 ym in fields with temperature gradients 
of the order of 5K/mm (cf. Walker et al., 1979) A common exam- 
ple of the phenomenon is the blackening of the glass globe 
of a kerosene lantern; the temperature gradient established 
between the flame and the globe drives the carbon particles 
produced in the combustion process towards the globe, where 
they deposit Thermophoresis also causes the dirt pattern on 
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the ceiling of kitchens 

The phenomenon of thermophoresis is of technical importance 
as IS discussea in detail by Green and Lane (1964) and Fuchs 
(1964). There are numerous technological and environmental 
problems wnich involve deposition of aerosol particles. Of 
current concern is the deposition of radioactive aerosol parti- 
cles insioe a nuclear reactor containment following a hypo- 
thetical accident, this plays a critical role in reducing the 
amount of raoi cacti ^ri ty that might subsequently escape to the 
atmospnere 

Thermophoresis has also been cited as a cause for the depo- 
sition of particulate matter on heat exchange surfaces with 
the attendant reduction of the heat transfer co-ef ficient . 
Due to the sharply incrcasj.ng cost of high-grade fuels, there 
is renewed interest in the combustion of coal- derived or lower 
grade 'residual* petroleum - derived gas turbine fuels contain- 
ing inorganic impurities. At the same time, there is a need 
for higher turbine inlet temperature to raise thermodynamic 
efficiency. However, the metal temperature is limited by struc- 
tural considerations and is maintained by film or transpiration 
cooling etc As a result, the high temperature gradients, and 
the large concentration of contaminants (ash, salts) present 
in the combustion products of lower grade fuels, leads to the 
formation and possible deposition of particles on gas turbine 
blade surfaces giving rise to highly accelerated blade corrosion 
and /or fouling rate. Hence, blade surfaces lose aerodynamic 
efficiency even if the changes in shape are small. Under these 
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combined circumstances, thermophoresis plays an important role 
in determining mass deposition rates as demonstrated by Goren 
(1977) and by Gokoglu and Rosner (1986) Such results and tech- 
niques are also necessary to aevelop systems for "clean" hydro- 
carbon fuels combustion since 1) local fuel pyrolysis and inade- 
quate oxidation rates caf^ cause soot particles production and 
aeposition, and/or 2) the air ingested for combustion may itself 
certain inorganic substances (e g , dust), which can deposit 
on gas-turbine components. 

Thermophoresis is also of practical use in removing aerosol 
particles like dust, mist, smoke etc from air in a room or 

from gas streams in various channels The thermal precipitator 
is one such piece of equipment, designed to utilize this pheno- 
menon positively It is sometimes more effective than electro- 

static precipitators in removing submicron - sized particles 
from gas streams (Talbot et al. , 1980). It is particularly 

effective in removing small particles with a diameter a micron 
or less Other important application of thermophoresis are 
in determining exhaust gas trajectories from combustion devices, 
in the manufacture of fumed silica, carbon black, and titania 
particles for paint industry, in the design of aerosol sampling 
devices, in chemical coating of metals, etc It has also been 
demonstrated recently that thermophoresis is the dominant mass 

transfer mechanism in the Modified Chemical Vapour Deposition 
(MCVD) and Outside Vapour Deposition (OVD) processes employed 
in the fabrication of optical fiber preforms (Simpkins et al. , 
1979; Walker et al , 1980, Alam and Mehrotra, 1987) It is 



by now apparent that silicon based optical fibers having low 
loss and wide band width are in the process of revolutionizing 
the technology of communications by replacing conventional 
copper cables 

The fact that scrubbing in aerosol scrubbers is more effec- 
tive when the dusty air is preheated has been ascribed to thermo- 
phcresis The reverse effect, namely, the repulsion of particles 
from a hot wall and the formation of the "dust- free" regions 
in the vicinity of hot walls, has also been observed Because 
such "clean" layers shield hot surfaces from particles in 
tne environmenr , the phenomenon has recently been recognised 
as an efficient means to protect electronic components from 
particle deposition during their manufacturing. 

1 1 Literature Review 


The basic physical processes involved in the phenomenon 
of thermophoresis were first investigated by Maxwell (cf. Ken- 
nard,1938), in an attempt to explain the radio-metric effect. 
Subsequently, this phenomenon has been the subject of theo- 
retical and experimental studies by a variety of investigators. 
The extensive literature on these investigations contains a 
number of conflicting results. A brief review of the existing 
important literature follows below : 

1.1.1 Thermophoretic Coefficient 


Small particles suspended in gas streams containing a tempe- 
rature gradient are known to be subjected to the thermophoretic 
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force and move steadily from hotter to colder regions in balance 
with VISCOUS drag As a consequence, the velocity acquired 
by the particles relative to the gas is known as the thermo- 
phoretic velocity Theoretical analyses yield the following 
linear relation between the thermophoretic velocity, w^ of 
a particle and the local temperature gradient, V t 

w^ = - K Vt ( T . T ) 

where t is amoient aosolute temperature of the gas, v is the 
kinematic viscosity of the gas, and K the dimensionless thermo- 
Dhoretic coefficient This coefficient is an important parameter 
in the particle transport and deposition The value of the 
thermophoretic coefficient depends mainly on the regime of 
flow as well as on the ratio of thermal conductivities of the 
particle and gas For free molecular flow 1), Waldmann 

found the coefficient to be 6/8 when there are elastic colli- 
sions of the molecules with the aerosol particles, and 6/(8 
+ IT ) for collisions resulting in condensation and reevaporation 
(cf Goren, 1977) For the slip flow regime (K^<<1) Epstein 
was first to calculate the thermophoretic velocity by solving 
the equation of viscous hydrodynamics (cf. Derjaguin et al , 
1976 ). With an arbitrary value of thermal slip coefficient, 
C = 0 . 75 , the Epstein's formula for thermophoretic coefficient 

O 

IS 
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where 

and k 
g P 

are the thermal conductivities 

of 

the 

gas and 

the 

particles , 

respectively Brock 

(1962) made 

the 

f ormula 

more 

precise to 

the technologically 

important 

transition range 

where 

the mean 

free path, x is of 

the order 

of 

the 

particle 

radiu 

s, a His formula for K is 






^ l+Ct X/a kp/kg 

^ (1 + i:p/2ig + Cp A/a kp/kg) ( ^ + 2C^ A/a) 

(1 3) 

where Cp, are temperature jump coefficient and velocity 

Jump coefficient, respectively and have values Cp = 2 16, 

= I 23 Broct's result, with = 0 75, is found to be not 

jin gooa agreement with experiment for particles of high thermal 
conductivity, although the discrepancy is much less than that 
obtained using the Epstein result 

An expression for calculating the thermophoretic coefficient 
has been proposed by Derjaguin et al (1976) The thermophoretic 
velocities measured by them support the expression 


(1 


A/a 


K 


.2 S- 


(1 +y^/2 ^ + C, . A7a k^/k„) 

P g t p g 


(1 4) 


where C = 1 1, close to theoretical value and C,= 2 17- Equation 
(1.4) IS the formula most frequently cited in the current lite- 
rature But because of the questionable nature of the theory 
underlying this result, plus the adjustment of the constant, 
it can hardly be regarded as more than an empirical formula. 


More recently, a critical review of the data and theore- 
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tical expressions for thermophoretic coefficient has been given 
by Talbot et al. ( 1980) They give evidence in support of their 
belief that the following expression for K is tolerably accurate 
for all regimes from free molecular flow to continuum flow 

-C, a/X 

(k /k + X/a) [1 + X/a (A,+ B,e ' ) ] 

1/ _ p p s B t ^ 1 ^ ^ ^ 

( 1 + 3 c X/a) (1 + 2 k /k + 2 C, X/a) 

m S P ^ 

(1.5) 

in which the constants have the values = 1 20, = 0.4l, 

= 0 86, Cg = 1 17, = 2 18 and C^= 1 14 The semiempi- 

rical relation (1 5) appears reasonably satisfactory for the 
entire range, 0 4 . However, accurate measurement of K 

IS difficult and the agreement between the available data and 
the theoretical values is not close Nevertheless, in most 
of the literature reviewed in section 112, it was contended 
that K can be treated as a specified constant of order 1 and 
independent of temperature in the flow field. 

1.1 2 Thermophoretic Flow Investigations 

Thermophoretic deposition of aerosol particles in flow 
fields have been the subject of many investigations. Laminar 
thermophoretic flow over a flat plate at zero incidence was 
analysed theoretically in detail by Goren (1977) for both cold 
and hot walls. Compressibility effects were retained because 
thermophoresis frequently occurs under large temperature diffe- 
rences for which the fluid properties may vary. The governing 
boundary layers equations were transformed to forms suitable 
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for solution by invoking Howarth -Dondnitsyn transformation. 
The similarity solutions with and without viscous dissipation 
were obtained in the absence of Brownian diffusion whose effect 
IS only to form a very thin sublayer even on the boundary layer 

scale For cold walls, the particle concentration drops monoto- 

nically as the wall is approached and a finite concentration 

at the wall is found depending upon the wall temperature For 
hot walls, the presence of a critical layer within the concentra- 
tion boundary layer is indicated by Goren He also emphasized 
the importance of the factor P^K in determining the concen- 
tration profile over a hot wall, where is the Prandtl number 
The effect of viscous dissipation on particle concentration 

profile, exhibited for cold wall conditions , indicates that 
wall particle concentration reduces with increasing dissipation 
rates 

Thermophoretic deposition of small particles in laminar 
incompressible flow through a tube with a prescribed wall tempe- 
rature distribution has been studied theoretically by Walker 
et al. ( 1979 ) who included the effects of Brownian diffusion 
by the use of matched asymptotic expansions They found that 
only a fraction of the particles present initially in the stream 
get deposited on the cold tube wall. This was because of the 
ultimate relaxation of the temperature gradient far downstream. 
It is evident that their similarity solution is valid only 
within a short axial distance from the start of the constant 
wall temperature. This solution was extended downstream 

using numerical methods It was also observed that cumulative 
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efficiency rapidly approaches an asymptotic limit for any wall 
temperature and value (along the axial direction) 

An analytical study of thermophoretic deposition of small 
particles was carried out by Homsy et al. (198I) This was 
the first investigation which dealt with an external flow in 
the presence of a non-zero pressure gradient. Their analysis 
was based on the development of thermophoretic flux in terms 
of Blasius series The values of particle flux coefficient 
were calculated numerically as functions of K and wall tempera- 
ture for fixed Prandtl number, Pj^= 0 71 Knowing the outer 
potential flow constants, one can obtain the corresponding 
deposition profile from these particle flux coefficients The 
local depositions profile plotted for transverse flow past 
a circular cylinder as an illustrative example, exhibits maximum 
deposition rate of the stagnation point and zero near the cylin- 
der angle, <i> = 90° For still larger 4) physically impossible 
negative values result because the truncated Blasius series 
IS not accurate The total deposition rate was found to decrease 
monotonically with increasing wall temperature. 

Thermophoretic transport of small particles through both 
laminar and turbulent natural convection boundary layers adja- 
cent to cold vertical plates have been analysed by 
Epstein et al. (1985). Introducing conventional Pohlhausen 
similarity transformation, the governing boundary layer equa- 
tions are reduced to ordinary differential equations by exploit- 
ing the concept of local similarity. Numerical solutions of 
these equations have been obtained with K and tw/t^ as 
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parameters Dimensionless particle concentrations at the wall 
as a function of K, for = 0 73, have been plotted for 

laminar flow It is evident from this graph that the wall parti- 
cle concentration decreases with decreasing tw/t^ and increasing 
K Interestingly enough, it is noted that while the particle 
concentration distribution for laminar flow differs distinguish- 
able in shape from that for turbulent flow, at the wall two 
profiles converges to practically the same particale concentra- 
tion 

/i detailed theoretical analysis of thermophoretic deposi- 
tion of particles on a cold surface has been presented by Bat- 


chelor and Shen (1985) 

By 

neglecting viscous dissipation and 

Brownian 

diffusion, a 

few 

exact relations 

are derived 

which 

are of 

application 

to different 

flow fields 

These 


relations suggest that the dimensionless wall particle concentra- 
tion Nw/N^ is approximately uniform over a cold isothermal 
Doundary (at wall temperature, tw) , that the total particle 

flux to the boundary is proportional to the total heat flux 
to the boundary, and that the relation between Nw/N and tw/t 
IS approximately the same for all flow systems These results 
are further compared with the available numerical results for 
a number of steady flow systems and found to be in good agree- 

ment except when both tw/t ^ and K are small 

Thermophoresis of aerosol particles in two-dimensional 

laminar flow in the presence of non-zero pressure gradient 

has been analysed in detail by Jayaraj (1988) The boundary 
layer-type governing equations are solved using a simple finite- 



difference marching technique for cold, hot and adiabatic wall 

conditions. While a non-uniform grid is used in the marching 
direction, the grid size in the transverse direction is self- 
adapted Results are presented for flow over an inclined plate 
and a circular cylinder for different values of tw/t^ keeping 
and ¥ as constant Irrespective of the geometry of flow, 

the particle concentration at the cold wall is found to be 
almost, independent of location and also of the angle of incli- 
nation (of the inclined plate). For both the inclined plate 
and tne circular cylinder, the thickness of the critical layer 
increases as one moves away from the stagnation point The 

influence of viscous dissipation on the particle concentration 
IS also discussed 

Quite recently, the thermophoretic deposition of small 
particles on cold surfaces has been treated by Shen (1989)- 
He investigates the external flow problem and proposes a general 
method of calculating the thermophoretic deposition rate of 
particles suspended in two-dimensional or axisymmetric flow 

past a body which is held at a lower temperature than the flow. 
Special consideration is given to the flow near the stagnation- 
points of a cylinder and a sphere, and the flow past a flat 

plate to show explicitly how the wall particle concentration, 

c /c and deposition rate are influenced by the density varia- 

w 0 ° 

tions. A general numerical procedure for this problem is used 
for calculations of the thermophoretic deposition of particles 
in flows past a cold cylinder and a cold sphere. Results for 
c /c at the stagnation point of the cylinder indicate that 

W 00 
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this quantity de^TPases with decreasing and increasing 

P K It IS also seen, for a range of values of P K and t /too , 

that c^/c^ remains almost constant along the surface This 
last result is also true for the flow over a sphere 

1.2 Objective of Present Investigation 

Very few investigators to date have considered the effect 
of gas thermophy sical property variations in the analysis of 
thermophoresis Only Goren ( 1977) and Shen ( 1989) took into 
accou’^t tne effect of density variation on the thermophoretic 
deposition Goren treated the thermophoresis of aerosol parti- 
cles in the laminar, compressible, boundary layer flow past 
a flat plate by using similarity transformations However, 
most of the practical applications of thermophoresis in exter- 
nal flow involve a non-zero pressure gradient (either due to 
the shape of the body or due to its orientation with respect 
to the main flow) for which no similarity solution is possible 

Shen studied thermophoretic deposition of small particles 
onto cold surfaces of bodies in two-dimensional and axisymmetric 
flows He included the compressibility effects with the assump- 
tion that the product yp is a constant where, y is the dynamic 
viscosity and p , the density of the gas. (This assumption 
was made by Goren as well) This is a well-known simplification 
and IS close to reality for not-too-large variations in tempera- 
ture. However, it should be emphasized here that thermophoresis 
IS frequently found under conditions of large temperature gra- 
dients. Thus, the application of above assumption may lead 
to erroneous results. In addition to this assumption, Shen 
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neglects the effect of viscous dissipation by taking M^<<1 which 
excludes most of the thermophoretic flows of engineering inte- 
rest For example, current combustion turbine blades experience 
high subsonic Mach numbers say, M^= 0 8 (cf Gokoglu and Ros- 
ner, 1986) It is v.ell-known that viscous dissipation modifies 
temperature profiles in momentum diffusion boundary layer 
Therefore, the assumption of no viscous dissipation will give 
inaccurate results There have been relatively few investiga- 
tions IP which the effect of viscous dissipation is included 
Moreover, Shen analysed the flow past circular cylinder by 
assuming inviscid, incom*pressible flow outside the boundary 
layei^ which differs very much from the real outer flow except 
for low Mach numbers 

The objective of the present work is to study the effect 
of variable thermophysical properties on velocity, temperature 
and concentration profiles with special reference to thermopho- 
retic deposition Accordingly, a pressure distribution is impo- 
sed on the laminar compressible boundary layer over a cold 
cylinder The pressure gradient is known from inviscid flow 
solution in the free stream The governing partial differential 
equations are solved by an implicit finite difference marching 
technique; while prespecified variable mesh size is used in 
the marching direction, the grid size in the transverse direc- 
tion IS self-adapted. The numerical solution includes the effect 
of viscous dissipation and neglects Brownian - diffusion (which 


influences only an extremely thin layer near the wall) 



CHAPTER 2 


MATHEMATICAL FORMULATION 

The first step in assessing the quantitative importance 

of t bermophoretic deposition of aerosol particles is to develop 
a mathematical model of the flow systems under consideration 

The actual flow fields encountered in most of the thermopho- 
retic applications are of a very complex nature The non-zero 
pressure gradient in the boundary layer over a body, common 
ir these applications, is due to either the shape of the body 
or to its orientation with respect to the undisturbed flow 
This pressure gradient causes further complications Certain 

idealizations of the physical model need to be made before 
a mathematical model can be obtained In this chapter, the 
governing equations fc- thermophoretic flow past a cold circular 
cylinder will be analysed 

2 i Physical Model 

Consider a stream of hot gas, containing suspended aerosol 

particles, flowing past a circular cylinder maintained at a 


lower 

temperature than that of the 

oncoming stream 

[ see 

Figure 

(2.1)] 

The temperature 

gradient 

established in 

the 

forward 

moving 

f orced-convection 

boundary 

layer between 

the 

surf ace 

of the 

cylinder and the 

gas-particle mixture drives the 

parti- 


cles towards the surface, where they deposit. Here, we consider 
the case of two-dimensional subsonic flow with significant 
Mach number effects. 


The particle concentration is assumed to be sufficiently 
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dilute so that the velocity and temperature profiles of the 
gas can be determined as in a flow free of particles and the 

thermophysical properties of the mixture are those of the gas- 
phase For particles of unit density, (i e., particle mass 

density comparable to that of the gas) and 1 y m radius, this 
assumptions limits the analysis to aerosol concentration less 

7 

than about 10 particles per cubic centimeter of the gas (Goren, 
1 97 7 ) This IS indeed the case for roost of the thermophoretic 

applications Additional assumptions usually made for the idea- 
lizatio! of a physical model in thermophoretic analysis (Goren, 

1 977, Walker et al , 1979; Epstein et al , 1985) are : 

(i) The particles are small enough so that in the absence 
of thermophoresis, the particle moves with the local 
gas velocity 

(ii) Thermophoretic velocity component in the streamwise 
airection is negligible 

(ill) The effect of Brownian diffusion is restricted only 
to an extremely thin layer next to the surface 

(iv) The effect of thermal radiation is negligible (This 

excludes the case of absorbing aerosols 

In addition, the thermophysical properties of the gas are 
assumed to be functions of temperature alone. 

2.2 Governing Equations 

The boundary layer equations governing the velocity and 
temperature fields have been thoroughly reviewed (White, 197^, 
Schreier, 1982; Schlichting, 1979). For a steady, laminar. 
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two-dimensional, compressible flow with variable thermophysical 
properties, these equations are (Schl ichting , 1979) • 


Total mass 
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where, (u,v) are velocity components in (x,y) directions [see 
Figure (21)], p, p and t are, respectively, the density, 
pressure and temperature of the gas, and p , C and k are 
the dynamic viscosity, specific heat at constant pressure, 
and thermal conductivity of the gas 

In the equations (2.2) and (2 3) the local pressure field 
IS assumed to vary only in the streamwise direction, i e. , 
pressure remains a constant across the boundary layer and is 
given by the inviscid flow solution More details will be given 
in Section 2.3 

The equation system (21) to (2.3) is mixed parabolic/hyper- 
bolic and, consequently, requires the following boundary condi- 


tions . 
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u (x , 0) 
V (x , 0) 
u (x , CO ) 

u (0, y) 
t ( X , 0 ) 

t (x , <») 
t (0, y) 


= U^Cx) . 

= u^(0) , 

= t, for isothermal wall, 

w 


(2.4) 


wnere is the free stream temperature, t^^, the surface tempe- 
rature and u (x) is inviscia flow velocity at the body surface 

It should be pointed out that the boundary layer equations 
are really not valid at x = 0 In computation, the effect of 

the singularity at the stagnation point is handled by selecting 
a very small step size in the x-direction (Hornbeck, 1973) 


Due to the boundary layer behaviour, the temperature gra- 

dient, — 2 — in the direction normal to the surface is very 

9t 

much larger than in the tangential direction, — . Also, 
the gas velocity component, u in the tangential direction is 
large compared to v, the gas velocity component in the normal 
direction This implies that the normal thermophoretic velocity 
components, v^, only needs to be considered. Using equation 
(1.1), v^ IS given by 


V 


t 



9y 


(2.5) 


As mentioned in section 1.1 1, Talbot et al. ( 1980) proposed 
a semiempirical formula for K which seems to be accurate for 
a wide range of Knudsen numbers. However, in most of the thermo- 
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phoretic calculation K is usually taken as a constant in the 
range of 0.2 to 12 For the present study, K is assumed 
as specified constant 

For compressible flow it is more convenient to introduce 
c, the particle concentration per unit mass of gas The particle 
concentration per unit volume, N and c are related by 

c = N/p (2.6) 


where p is the density of the carrier gas 

With the usual boundary layer approximation, the species 
mass conservation equation for the physical model depicted 
in section 2.1 can be written as 

( PCu)+-^[ po(v+ v^)l = -j|-[ (27) 

where, D is the Brownian diffusivity of the particles. Equation 
(2.7) can be rewritten in simplified form using the total mass 
conservation equation (2 1) as 

pu-^ + PV-|§- +-3|- (pov^) = •^[ D -^(PO )] (28) 

The boundary conditions are 

c (0, y) = , 

c (x, °°) = , 

and c (x, 0) = 0 (2.9) 
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if the wall IS totally absorbing Here, c^is the particle 
concentration in the free stream 

The Brownian diffusivity, D, of aerosol particles of dia- 
meter 0 05-1 0 ym is very much smaller than any of the gas 
dif fusivities (for mass, momentum, heat). For example, for 
particles with a size of 0 1 ym the particle Brownian diffusi- 
vity IS of the order of lO'^cm^s””' Thus, for particle of this 
size or larger the Brownian diffusion sublayer is so thin even 
in comparison with the boundary layer thickness that it does 
not affect the thermophoretic deposition rate (Goren, 1977; 
Walker et al. , 1979) The only effect of Brownian diffusion 

IS to create a thin particle concentration boundary layer adja- 
cent to the cold surface. Hence, the particle diffusion term 
on the right- hand side of Equation (2 8) could be dropped 
and can be rewritten as 


P u 


3c 

3x 


P V 


3c 

3y 


3y 


( p cv. ) = 0 


( 2 . 10 ) 


Equation (2 10) is of the first order in y . It implies that 
only one y-boundary condition (at infinity) will be needed 
for its solution. In the absence of diffusion, a finite particle 
concentration at the wall is expected for the cold wall condi- 
tion. This wall particle concentration may be obtained by solv- 
ing Equation (2 10). In fact, this concentration corresponds 
to the outer edge of the very thin Brownian diffusion sublayer 
adjacent to the cold surface. 
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2 3 Pressure Distribution 

The calculation of the pressure gradient in the boundary 
layer is essential for the solution of the momentum and energy 

equations The pressure and velocity distribution at the wall 
are calculated on the basis of inviscid flow theory and then 
used as conditions at the outer edge of the boundary layer. 
Since the thickness of the boundary laj er is not known a prion, 
the conditions are applied at y = in the solution In short, 

the condition at y = 0 for inviscid flow are assumed to prevail 

at y = “in the viscous flow problem. This procedure is, of 
course, widely used in Fluid Mechanics, its justification being 
that it is found to give good results in thin boundary layers 
Solutions found in this way show that the velocity reaches 
99? of its inviscid value at a short distance from the wall, 

thus confirming to the boundary layer concept The pressure 
in the bounaary layer is found from (Schreier, 1982) 


1_ _d£ 

P d> 


u 


d u 
dx 


s 


( 2 . 11 ) 


The pressure gradient, — is usually computed by consider- 
ing inviscid flow over the body surface Two-dimensional, invis- 
cid, compressible, subsonic flow over a circular cylinder is 
assumed to be uniform at infinite with zero strength and zero 
circulation. The following expression for dimensionless velocity 
over the circular cylinder is obtained by the Rayleigh - Janzen 
method (Schreier, 'S982) 
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(2 12 ) 


where L is the radius of the cylinder and x corresponds to 
the location on the cylinder surface from the stagnation point 
The equation (2.11) can be rewritten with the help of equation 
(2.12) as 
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[2 + ^ ) - 6 cos^ ( x/L) 


(2.13) 


2 4 Constitutive Relations 

There are only six equations in nine unknown variables, that is 

u, V, V. , p, t, c, p, p, and k In order to solve these equa- 
^ S 

tions, three more constitutive equations of state will be nece- 
ssary. One such well-known constitutive equation is equation 
of state for perfect gas, i.e., 

p = P R t (2 14) 


The dynamic viscosity, y, and thermal conductivity, k of the 

o 

gas are generally assumed to be functions of temperature only. 
The commonly employed power-law relationship will be used here. 
These expressions for y and k are (cf. Gokoglu and Rosner, 

o 

1986 ) : 


y/y 


t ^0.70 
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(2 15) 
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(2 16) 


wnere u and k are the values at tne stream temperature t 

^00 S °° QO 

2 5 Di’^ensionless Representataon 

Thp dimensionless form of system of equation are more compact and have 
les. Steen velocity, temperat’ure and concentration gradients The sharp 
gradients usually create convergence problems in numerical 
calcu j ati O’^s In addition, nondimensionalising of the governing 
ecuat loni- leads to the identification of nondimensional numoers 
bas'^d or- wn^ch tne result can be preaicted Following dimension- 
less ■variaole are introduced 
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Here, L the radaus of the circular cylinder and R is the 

Gqo 

frte ctreani Re/rolas number defined as 
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( 2 . 18 ) 


III additic’', the following nondimensional numbers are introduced 
in equation • 


: f O’' cti ’"'uirber , P 
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(2 19) 


wnere, y the specific heat ratio and R, the gas constant. 

l£..rig abov-e dimensionless variables and nondimensional 
riumbcrs, the system of Equations (2 1) to (2.5) and (2 9) to 
(2 16) can oe rewritte’^ in dimensionless form as given below 

Total macs eonscr \/at ion : 
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( 2 . 24 ) 


(2 25 ) 


( 2 . 26 ) 
(2 27 ) 
(2 28 ) 


and 


( 2 . 29 ) 
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? 0 otratfaKv for Solution 

governing equations (2 1) to (2.3) and (2 10) are 
n'lnear coupieci partial differential equations They are 
sj Cv^mol ^eatea that an analytical solution especially in the 
r'r’f-^scnee of nor-zero pressure graaient is not only difficult 
-ot prooacly n^nossible Hence it is essential to seek a numeri- 
cal s. -utioi lo these ecuation 

The solution to these system of equation can be obtained 
jC.rg the f 1 r I r e-dif ference technique while marching in the 
y-iz" ‘tion. Tne oressure gradient required for the solution 
" mo^'untum arc energy equation can be computed by inviscid 
cumpressible flow considerations. Both streamwise momentum 
ctj'd total mass conservation equations can be solved simultane- 
ously in an iterative manner to arrive at the correct velocity 
profiles Once the two velocity components are known, the tempe- 
rature profile can be obtained by solving the energy equation 
Finally, the particle concentration equation can be solved 
using the known values of density and dynamic viscosity at 
the current location. 



CHAPTER 3 


COMPUTATIONAL TECHNIQUE 

Tfi'? prjf'bs uf obtaining the computational solution con- 
5ist,c of tw.' stages The first stage converts the continuous 
pa’^Mil 01 f fo*’t.r'tial equations and auxiliary (boundary and 
i^it^al' oc^'dat^ors into a discrete system of algebraic equa- 
tionr T‘ 1 = stage is called ’discretisation’ The second stage 
»^equi’'o: a.r- 'tcuatior sol/er’ to provide the solution of the 
s,'''tt”‘ of algebraic cauations 

’ umber of cncices are available for the discretisation 
of I ru g'V'j>'‘ning partial differential equations • the finite- 
aifference method, the finite-element method and the spectral 
method The f ini te-dif ference method is used here The finite- 
diffe^’ence methoa consists of the construction of a discrete 
grid, the replacement of the continuous derivatives in the 
gr\tr-nirjg pamal differential equations with equivalent finite- 
differe’^ce expressions and the solution of the resulting algeb- 
raic equations by an algorithm This chapter is concerned with 
the computational techniques for obtaining and solving the 
system of algebraic equations generated from the governing 
equations (2 20) to (2.25). 

3 1 Discretisation of the Governing Equations 

In local methods, like the finite-difference method, 
the discretised algebraic equations link together values of 
the dependent variables at adjacent grid points In order to 
convert the dimensionless differential equation (2.20) to (2.25) 
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j li f J M t f ^.ffertr'ce form, a two-dimensional rectangular mesh 
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art.' gr.jd irit-r’vals That is, (0,0) represents the stagnation 
f t 1 ” terms of spatial coordinates and finite-difference 

-r dices The cylinder surface is represented by k = 0 and the 
‘^s.re cf toundary layer by k = n+ 1 

P'Cc'U-.t" tr.p governing equations are of mixed hyperbolic/ 
p --’’'c'b'hi ic tvp<', The comoutat ion can be carried out by "marching” 
fr >w^.fd -r, tv e /-direction The momentum, energy and concentra- 
tif-n equations are solved implicitly at every step The total 
mctSsS oonser Vet 1 on equation is solved simultaneously with the 
momentum equation by an iterative algorithm Forward differenc- 
ing in tnc X-direction and the central differencing in the 
Y-dirertion are used in Equations (2 21) and (2 22) Forward 
differencing is employed in Equations (2.20) and (2.24) The 
f inite-aif f erence form of these equations, in a form readily 
adaptable for solution, is given below . 


Total mass conservation equation 
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The rmopho retie Velocity : 
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( 3 .^) 


Constitutive Relations : 
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^ Vl,k 

j-1.0 

(3.5) 
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The above discretisation is 

second order accurate in 

Y for 


U and T. 

In order to satisfy the sufficient conditions 
for a meaningful solution of a triadiagonal set of equations 
(Roache, 1982) different discretisation forms are used for 
species mass conservation equations (2.25) depending upon the 
signs of V . These equations are given below : 


Particle Concentration Equation : 
Cold wall ( V > 0 ) ; 
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Cold Wall ( V < 0 ) 
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This discretisation of species mass conservation equation is 
first order accurate in AX and AY. It should be emphasized 
here that V does not necessary change sign at all X-locations. 

3.2 Method of Solution 

The computation of the values of U, V, T and C is carr- 
ied by "marching”, starting from the stagnation point (X=0). 
The equations (3.1) to (3-3) are linearised to facilitate the 
search for the solution. This is done by taking the values 
of p , p , and k (as the case may be) at the location (j+1) 

O 

equal to the corresponding values at the location (j). In addi- 
tion, the coefficients ^j+1,k non-linear 

terms) in Equation (3.2) is replaced by their computed values 
in the previous iteration, indicated by the primed quantities 
in Equation (3.2). The linearised tridiagonal set of momentum 
equation is solved iteratively using the discretised continuity 
equation (3.1) to update V at the location (j+1). 

Once the values of U and V have been computed at the 
(j+1) location, the discretised energy equation (3-3) written 
for k = 1(1 )n becomes a tridiagonal set of h linear equations 
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that can be easily solved for n unknown values of T at the 

location (j+1). 

The temperature values computed at the (j+1) location 
are substituted in constitutive equations (3-5) to (3.7) to 
obtain n values of p*, y*, and k* at (j+1) location. The n 

O 

unknown values of V,p at this location (j+1) can now be obtained 

using the recent computed values of p , u , and k at the loca- 

S 

tion (j+1) using equation (3.4). 

The solution of the discretised equation (3.8) or (3.9) 
(depending whether V is positive or negative) can be obtained 
by using the computed values of U, V, p* and V,j, at the (j+1) 

location. In the case of positive V, Equation (3.8) is selected 
as the discretised concentration equation. However, this equa- 
tion can not be solved at ksO, in which case Equation (3-9) 

( multiplied by a minus sign to preserve diagonal positiveness) 
has to be used for this point. This procedure yields a tridia- 
gonal set of (n+1) linear equations for the (n+1) unknown values 
of C at the (j+1) location. In the case of negative V, Equation 
(3*9) written for k=0(1)n leads to a bidiagonal set of (n+1) 
linear equations for the (n+1) unknown values of C. 

3.3 Computational Details 

The velocity, temperature and concentration boundary 
layer thicknesses increase as we move along the surface of 
the cylinder away from the stagnation point. Thus, the three 
boundary layer widths are adjusted separately, in the compu- 
tation, while marching in the X-direction. This is done by 
increasing the number of mesh points, n, in the Y-direction, 
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at each X-location. The value of n is chosen so as to ensure 

that there are at least 3 grid points where U kU , T k 1, and 

s 

C . This is to ensure that the computation indeed spans all 
three boundary layers. 

The step size X is increased systematically after march- 
ing certain number of steps in the X-direction. This has been 

done to improve computational efficiency. Very small tx (of 

-5 

the order of 10 ) is chosen while starting from the stagnation 

point. This value is almost doubled at an interval of about 
10 marching steps. However, beyond X=0.1, the same AX is used 
for increasingly large number of steps. Thus, in the span of 
X from 0 (stagnation point) to 1.60 (far away from the stagna- 
tion point), twelve different step sizes are used for computa- 
tion. Table 3*1 gives the value of AX used and the number of 
steps marched with each AX. 

Because central differencing is employed in the Y- direc- 
tion, it is preferable to keep AY uniform to preserve second 
order accuracy. But then the number of simultaneous equations 
to be solved increases excessively because the steep gradients 
near the wall restrict AY to a very small value. This increases 
computer time, storage space and also round-off error. A practi- 
cal solution to this problem is to employ a fine mesh near 
the cylinder surface and a relatively coarse mesh away from 
it. Accordingly, a self- adaptive grid scheme is used along 
the Y-direction. This needs modification of the discretised 
equations ( see Hornbeck, 1973;Appendix) at locations where the 
mesh size in the Y-direction changes. The minimum value of ^Y 
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Table 3"1 : Marching Pattern along X-direction 


Serial No. 

ax 

Number of 
steps 
marched 

Final 

X value 

1 . 

0.00001 

20 

0.0002 

2. 

0.00002 

10 

0.0004 

3. 

0.00005 

12 

0.0010 


0.00010 

10 

0.0020 

5. 

0.00020 

10 

0.0040 

6. 

0.00050 

12 

0.0100 

7. 

0.00100 

10 

0.0200 

8. 

0.00200 

10 

0.0400 

9. 

0.00500 

12 

0.1000 

10. 

0.01000 

25 

0.3500 

11 . 

0.01000 

105 

1.4000 

12. 

0.01000 

20 

1.6000 
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IS set at 0.005 and the maximum at 0 05. Self-adaptive grid 
generation technique (see section 3-^) distributes values of A Y 
withan this range. This range of^AY was chosen, after numerical 
experimentation, to optimise the calculations with respect 
to computational speed and accuracy 

3.^ Self-Adaptive Grid Scheme 

A number of grid generation techniques are in use for 
the solution of the boundary layer equations. The aims of an 
adaptive grid scheme is to optimise the distribution of grid 
points while controlling grid skewness. The grid must remain 
as orthogonal as possible since excessive skewness causes incre- 
ased truncation errors. Most of the adaptive grid methods essen- 
tially attempt to equidistribute some positive weight function 5 
w(Y), of the solution over the field. That is. 


Y 

/ 

Y 


k +1 

k 


w(Y) dY 


constant 


(3 10 ) 


for an one-dimensional problem. In discrete form. Equation 
( 3 . 10 ) may be written as 

A Y. w r constant, (3-11) 

K K 

where AYj^ is the grid interval, and is the weight 

function (assumed to be constant over the interval). 

The self-adaptive grid scheme employed for the present 
computation is essentially the same as that suggested by Naka- 
hashi and Deiwert (1986,1987). This method redistributes the 



35 


is set at 0.005 and the maximum at 0.05. Self-adaptive grid 
generation technique (see section 3-4) distributes values of A Y 
within this range. This range of^AY was chosen, after numerical 
experimentation, to optimise the calculations with respect 
to computational speed and accuracy 

3. '4 Self-Adaptive Grid Scheme 

A number of grid generation techniques are in use for 
the solution of the boundary layer equations The aims of an 
adaptive grid scheme is to optimise the distribution of grid 
points while controlling grid skewness. The grid must remain 
as orthogonal as possible since excessive skewness causes incre- 
ased truncation errors. Most of the adaptive grid methods essen- 
tially attempt to equidistribute some positive weight function , 
m(Y), of the solution over the field. That is, 

^k+1 

w(Y) dY = constant (3 10) 

^k 

for an one-dimensional problem. In discrete form. Equation 
(3.10) may be written as 

AYj^ = constant, (3-11) 

where AYj^ is the grid interval, and is the weight 

function (assumed to be constant over the interval). 

The self-adaptive grid scheme employed for the present 
computation is essentially the same as that suggested by Naka- 
hashi and Deiwert (1986,1987). This method redistributes the 
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pt id optimally, reducing the overall solution error 

1 r: li'ldit inn it controls the orthogonality and smoothness of 
tJic gr id points. A deterministic expression for the weight 
f.-ii t i on has been suggested which provides minimum empiri- 

>. ar well as a good control of the parameters to the user. 

A' : ug.r-^sted by Nakahashi and Deiwert ( 1987), the expression 
f r w- igl/ function, chosen to be 

0)^, = 1 + A , (3.12) 


Wlj. 
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grid 
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r k min , . 

■‘K~f-f " tO-IOl 

max min 

tM* censtant A is zero, the grid points will be equally 
d The quantity fj^ is the non-negative driving function 
on the solution of the problem. This fj^ is usually taken 
‘ the gradient of the flow solution over the n grid inter- 
The quantities fj^ and are assumed to be constant over 
grid interval. Where A and B are positive constants The 
spacing IS obtained by 


“k oil 

whprp f and f are the maximum and minimum values of f, 
wnere arm r 

respectively and H is the length of the line in physical space, 
and n is the total number of grid intervals. The constant A 
controls the ratio of the maximum to the minimum grid spacing . 
This ratio is initially specified. Thus the constant A in Equa- 
tion (3.12) is given by ^Yj^ax ^^min 


A 


AY 


max 


AY 


min 


(3.15) 
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The constant B in Equation (3 12) is chosen such that the mini- 
■UDi grid spacing Mini a is equal to the specified A 

and this can be achieved by solving the equation 

►’^"1 \ I = (3-16) 

The Newton - Raphson method can be used to solve Equation (3-16) 

0 

starting from a first guess B, by tne iteration relation 


,m+ 1 


O® * "D® 

B + A B , 


(3-17) 
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( 3 . 18 ) 
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and represents the guess of B. The quantity 
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3Mim 1 A Y. 
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(3-1^) as 
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be evaluated using Equations (3*12) and 


3 Min I A Yj^ I 
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n 


Sil ‘ 


IL w. 


) f® In f . , (3.19) 

J J 


where and AYj^ correspond to Min 1 A Yj^ 1 . The function fj^, 

is interpolated on the new grid resulting from Equation (3.14) 
with the new value of B, and the whole process is repeated 
till a reasonably low value ofl AB 1 is achieved. 

A FORTRAN subroutine based on the above method controls 
the grid self-adaptation in the Y-direction. This subroutine 
is called at prespecified X-locations. The function f^^ can 
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be 3ny of the Y~gradients of U, T or C. However, for the present 
computation, the Y-gradient of U is chosen as it generally 
has the steepest slope. The minimum value of aY is set at 0.005 
and the maximum at 0.05. 

3.5 Validation of Computer Code 

In order to validate the computer code used, the flow 
past a flat plate at zero incidence with zero pressure gradient 
and of constant properties was computed. The results obtained 
are compared with the Blasius solution for the flat plate. 
These comparisons are depicted in Figures (3-1) and (3*2). 
The.se two plots correspond to two different streamwise posi- 
tions. Even at X = 0.01 (i.e. very close to the stagnation 

point) the two results are close to each other. The result 
obtained by computation matches vetyweEwith the Blasius solution 
at X a 1.57. 











CHAPTER It 


RESULTS AND DISCUSSION 


7tiv conruier code developed and tested for the flat 
plate ha.'j t-eer; ured to calculate the thermophoretic flow past 
a cold circular cylinder. Using the theoretically obtained 
distrituticn given by Equation (2.13), results are presented 
for a range of values of Mach number, wall temperature and 
thernscphorei;!' coefficient. Special emphasis is given to the 
results obtained by the consideration of variable thermophysical 
properties. The velocity, temperature and particle concentration 
profiles in the boundary layer are shown in Figures 4.1 to 
4.11 for Fj. s 0.7 at various angular locations on the cylinder. 
A comparison between the distributions obtained with constant 
properties and variable properties is illustrated in every 
case. These profiles are shown versus the nondimensional normal 
distance Y, at different locations along the cylinder surface. 


4.1 Hydrodynamic and Thermal Boundary Layers 


The hydrodynamic and thermal boundary layers over a cir- 
cular cylinder have been investigated extensively by various 
researchers (cf. Schlichting, 1979 and White, 1974). These 
results for velocity and temperature boundary layers presented 
here for the sake of comparison and completeness. 


4.1.1 Longitudinal Velocity Profilfg 




The dimensionleaa . 
files for variable 
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thc^ oy]ind€‘r are shown in Figure 4.1 These profiles seem to 
be in agreement with profiles shown in literature. Velocity- 
distribution obtained for the constant properties case, at cj) = 
^ presented in Figure 4.2, along with the corresponding 
profile for the variable properties case at the same location. 
The two velocity profiles are considerably different from each 
other. Notably, the profile for the variable properties flow 
has a steeper gradient compared to that of constant properties 
profile. 

4 1.2 Transverse Velocity Profile 

Figure 4.3 shows the transverse velocity profile 
at different X-locations for the variable properties 
case (Tw = 0.25, M = 0-3)- It is noted that the transverse 

velocity component is negative over a considerable portion 
of the cylinder surface (upto ^ = 7(f). The effect of variable 
properties on the transverse velocity distribution has been 
illustrated in Figure 4.4, again for Tw = 0.25, = 0,3. It 

appears that the assumption of constant properties leads to 
significant differences in the calculation of V at Tw = 0.25. 


4.1.3 Normalised Temperature Profiles 

Figure 4.5 shows the normalised temperature distri- 
bution in the thermal boundary layer for the variable proper- 
ties case and an isothermal wall condition (Tw = 0.25). This 
plot indicates that the temperature gradient for <|> = 1° is 

greater than that at 4> ~ 90°, as expected. 
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ti.TT0CT OT Properties vQriQxio n on the Longitudinol 
Velocity Distribution for Flow past a Cold 
Circular Cylinder (0 = 60°) 






igurc 4*4 Effect of Properties Voriolion on the Tronsversc 

Velocity Distribution for Flow past a Circular 
Cylinder 
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The effect of properties variation on the normalised 

temperaturt: distribution is depicted in Figure 4.6. Once again 
the profiles drawn for the constant properties case are quite 

different from that of the variable properties case for Tw 
= 0.2b. 

The effect of Mach number on the normalised tempera- 
ture profile is illustratedin Figure 4.7. It is noted that for 

a small change in the value of , the temperature profile does 
not alter appreciably. However, when M,o is increased from 0.2 
to 0.5 the temperature profile is modified to a great extent. 
It may be noticed that the = 0.5 profile has a less steep 
gradient than the lower Mach number profiles. One could expect 
that at higher Mach numbers the increased viscous dissipation 
would lead to higher temperature gradients. However, Gokoglu 
and Rosner have shown that even at M«>= 0.6, viscous dissipation 
hardly plays any role for Tw<0.6 (cf. Jayaraj, 1988). Figure 
4.6 is plotted for Tw = 0.25; therefore viscous dissipation 
can be assumed to have negligible effect in this case. Thus 
the profiles shown in the figure are probably influenced by 
the thermophysical properties variation at the various Mach 
numbers. However, we have not investigated this conjuncture 

in greater depth. 

4 . 2 Concentration Bo undary Layer 

f 

\ ^ f 

4.2.1 Particle Concentrat ion Profile^ 

0 • ■ • 

Figure il.a rtoOT tM perWel* profiles 

at various F-iocatlons. »lti^ the prof «•« stallar. 
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Figure 4-8 Particle Concentration Profiles over a Cold 

Circular Cylinder at Different X-Locations 
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they are not identical, i.e., they do vary along the surface 
of the cylinder. However , the particle concentrations at the 
wall are essentially constant throughout the entire range of 
X. This remarkable fact has already been noticed for the case 
of constant property (Batchelor and Shen, 1985, Epstein et al, 
1985 and Jayaraj, 1988) and for variable density flows (Goren, 
1977 and Shen, 1989). 

Figure 4.9 demonstrates that the concentration profi- 
les change significantly upon making the constant property 
assumption. In these computations too, the wall concentration 
IS essentially constant and independent of X. However the cons- 
tant property assumption over-estimates the wall concentration 
by about 20 % for this case (Tw = 0.25, = 0.7, M „= 0.7, 

K = 0.75). 

Table 4.1 gives a comparison of values obtained 

in the present study with those obtained by Shen (1989). 

Figure 4.10 displays the particle concentration profiles 
at a single X-position ( (t) = 60°) for flows of Mach numbers 

0.2, 0.3 and 0.5. It is clear that the Mach number has a notice- 

able, but not pronounced, effect on the concentration profile 

Figure 4.11 demonstrates the effect of wall tempera- 
ture on the particle concentration profile at a particular 
X-location ( <j> = 60°). The figure shows that the profiles become 
less steep, and the particle concentration at the wall increa- 
ses, with increasing wall temperature. This is similar to the 
results already obtained for the constant property case by 



igure 4*9 Effect of Properties Variation on the Particle 
Concentration Profiles for Flow past a Cold 
Circular Cylinder 



Values of Wall Particle Concentration near the Stagnation 
Point on a Cold Circular Cylinder (P = 0.72) 





0 


Figure 4-11 Effect of Wall Temperature on the Concentration 

Profiles for Flow past a Cold Circular Cylinder 
{cp= 60) 
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Jayara j ( 1 9B8) . 

The effect of the product P^^.K on the particle concen- 
tration at the wall, for different wall temperatures, is shown 
in Figure 4.12. This figure shows the following trends : (a) 

when P^.K increases, for a given wall temperature, the wall 
particle concentration falls and (b) for increase in wall tempe- 
rature, for a given ?^.K, the wall concentration increases. 
When P^.K is zero the wall concentration is the same as the 
particle concentration in the free stream. These results are 
similar qualitatively to those obtained by Goren (1977) for 
forced convection over flat plate (variable density), by Epstein 
et al. (1985) for natural convection over a vertical flat plate 
for both laminar and turbulent flow (constant property), and 
by Shen (1989) for forced convection over a cold cylinder (vari- 
able density). Thus the relationship depicted in Figure 4.12 
IS remarkably universal. 










CHAPTER 5 


CONCLUSIONS 

Thermophoretic effects in most industrial application 
occur in conjunction with high temperature gradients and signifi- 
cant (high subsonic) Mach numbers. Although this fact is well- 
known, little work has been published in literature, so far, 
which account for the variable properties in such flows. Many 
researchers have concentrated on the constant properties case 
and low Mach numbers. Those who considered high Mach number 
flows have included only simple formulations (and strong assump- 
tions) in their handling of the variable properties inherent 
in these flows. 

This thesis has considered the variable properties 
flow of significant Mach number over a cold cylinder. Although 
the results have a qualitative similarity with the previous 
studies, this numerical work for a standard flow has demonstra- 
ted that the variable properties, high Mach number case is 
quite diferent from the constant properties, low Mach number 
one. In particular, the velocity, temperature and particle 
concentration profiles show significant differences in the 
two cases. 

5.1 Recommendations for Further Work 

The field is wide open for studies incorporating 
other prooesees whioh also occur In situations involving thermo- 
phoresis. For example, bouyanoy effects which are particularly 
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strong in these situations, due to high temperature gradients, 
have not yet been studied in flows of engineering applications. 
Turbulence effects on thermophoresis can also be incorporated. 
No researcher has yet published any work that incorporates 
three-dimensional effects 



LIST OF REFERENCES 


1. Alam» M.K. and Mehrot^'a, S., 1987, ’Thermophoretic Depo- 

sition of Particles in Optical Fiber Preform Fabrication, 
ASME paper No. 87-HT-6, 2Hth Nat. Heat. Transfer Conf . . 

Pittsburgh, Aug. 1987. — — ■ 

2. Batchelor, G.K. and Shen, C. , 1985, 'Thermophoretic 

Deposition of Particles in Gas Flowing Over Cold Sur- 

Colloid Interface Sci. . Vol. 107, pp. 21-37. 

3 . Brock, J.R., 1962 , 'On the Theory of Thermal Forces Acting 

on Aerosol Particles', J. Colloid Sci. , Vol. 17. dd. 768- 

780 . 

M . Derjaguin, B.V., Rabinovich, YA.I., Storozhilova, A. I 

and Shcherbina, G.I., 1976, 'Measurement of the Coeffi- 

cient of Thermal Slip of Gases and the Thermophoresis 
Velocity of Large-Size Aerosol Particles', J. Colloid 
Interface Sci . . Vol. 57, pp. •451-'46l. 

5 . Epstein, M- , Hauser, G.M. , Henry, R.E. , 1985, 'Thermo- 

phoretic Deposition of Particles in Natural Convection 
Flow From a Vertical Plate' , ASME J. Heat Transfer , Vol. 
107, pp. 272-276. 

6- Fuchs, N.A., 1964, The Mechanics of Aerosols , Pergamon 

Press, New York. 

7 . Gokoglu, S.A. and Rosner, D.E., 1986, ' Thermophoretically 

Augmented Mass Transfer Rates to Solid Walls Across Lami- 
nar Boundary Layers', AIAA J . , Vol. 24, pp. 172-179- 

8. Goren, S.L., 1977, ’ Thermophoresis of Aerosol Particles 

in the Laminar Boundary Layer on a Flat Plate, 'J. Colloid 
Interface Sci. , Vol. 61, pp. 77-85. 

9 . Green, H.L. and Lane, W.R., 1964, Particulate Clouds : 

Dusts, Smokes and Mists , 2nd Ed., E. and F.N. Spon Ltd., 
London. 

10. Homsy, G.M., Geyling, F.T., and Walker, K.L., I 98 I, 'Bla- 

sius Series for Thermophoretic Deposition of Small Parti- 
cles', J. Colloid Interface Sci. , Vol. 83 , pp- 495-501. 

11. Hornbeck, R.W., 1973, Numerical Marching Techniques for 

Fluid Flows with Heat Transfer , SP-297, NASA, Washington, 
D.C. 

12. Jayaraj, S. , 1988, 'Thermophoresis of Aerosol Particles 

in Two-Dimensional Laminar Flow, Pb > 

Dept, of Mechanical Engineering, I.I.T. KANPUR. 

13 . Kennard, E.H. , 1938, Kinetic Theory of GM®®> McGraw- 

Hill, New York. 



62 


14 . 

15. 


Nakahashi, K. and Deiwert, G.S., 1986 

Adaptive Grid Method, * A1AA J. . Vol’ 


'Three Dimensional 
24, pp. 948-954. 


Nakahashi, K. and Deiwert, 
Grid Method with Application 
Vol. 25, pp. 513-520. 


G.S., 1987 , 'Self-Adaptive 
to Airfoil Flows', AIAA J. . 


16 . 

17. 

18 . 


Roache, P.J., 1982, Computational Fluid Dynarmns 

Ed., Hermosa Publishers 7 Albuquerque. ’ 


2nd 


Schlichting, H. , 1979, Boundary Layer Theory,, 7th Ed 

McGraw-Hill, New York. ^ 


Schreier, S. , 1982, Compressible Flow . John Wiley & Sons 

New York. ’ 


19 . Shen, C., 1989, 'Therraophoretic Deposition of Particles 

Onto Cold Surfaces of Bodies in Two-Dimensional and Axisy- 
mmetric Flows', J. Colloid Interface Sci.. Vol. 127. 
pp. 104-115. 

20. Simpkins, P.G. , Greenberg-Kosinski, S. and Macchesney, 

J.B. , 1979 , 'Thermophoresis : The Mass Transfer Mechanism 

in Modified Chemical Vapor Deposition', J. Appl. Phys., 
Vol. 50, pp. 5676 - 5681 . 

21. Talbot, L., Cheng, R.K., Schefer, R.W. and Willis, D.R., 
1980 , 'Therraophoresis of Particles in a Heated Boundary 
Layer’, J. Fluid Mech. . Vol. 101, pp. 737-758. 

22. Walker, K.L., Homsy, G.M. and Geyling, F.T , 1979, 'Ther- 

mophoretic Deposition of Small Particles in Laminar Tube 
Flow', J. Colloid Interface Sci ., Vol. 69, pp. 138-147. 

23 . Walker, K.L., Geyling, F.T., and Nagel, S.R., I 98 O', 

Thermophoretic Deposition of Small Particles in the Modi- 
fied Chemical Vapor Deposition (MCVD) Process', J. Amer. 
Ceram. Soc. . Vol. 63, pp. 552-558. 

24. White, F.M., 1974, Viscous Fluid Flow . McGraw-Hill, New 

York. 


AII0033 


T« 

T«3-e 


Dal* Slip 


Thif bool( is to b« rtturntd on th* 
date last stamped. 


I 



